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We present the results of the Model Unspecific Search in CMS (MUSiC), which systematically scans the data 
taken by the CMS detector for deviations from the Standard Model predictions. Due to the minimal theoretical 
bias this approach is sensitive to a variety of models for new physics. Events containing at least one electron 
or muon are classified according to their content of reconstructed objects (muons, electrons, photons, jets and 
missing transverse energy). A broad scan of three kinematic distributions in those classes is performed by 
identifying deviations from Standard Model expectations, accounting for systematic uncertainties. 
In this particular search data taken by CMS in the year 2010, corresponding to an integrated luminosity of 
36.1 pb _1 , have been analysed. 



1. Introduction 

The start-up of the LHC brings forth a new era of 
high energy particle physics. What we will find is yet 
unknown, however, there is a large number of theories 
predicting possible outcomes. Many of those theo- 
ries are tested by dedicated analyses at CMS and the 
other LHC experiments. However, new physics could 
as well manifest itself in ways no-one has thought of 
yet. For this purpose a Model Unspecific Search in 
CMS: "MUSiC" has been implemented. 

MUSiC automatically scans and compares the mea- 
sured data to the simulated SM expectations without 
assuming any specific model of new physics. Thus, 
covering a large phase space, it is sensitive to a lot of 
potential signals. 

Any uncovered significant deviation needs addi- 
tional interpretation, such that its origin can be de- 
termined. Possible causes could be insufficient under- 
standing of the collision, event generation or detector 
simulation, or indeed genuine new physics in real data. 
Thus, the output of MUSiC must be seen as only the 
first, but important step in the potential discovery of 
new physics. 

More details on this and following topics can be 
found in the P hysics Analysis Summary CMS Collab- 
oration 2011 . 



2. Implementation 

New physics usually shows up in distinctive final 
states and so the first task is to sort the events. We 
consider the following physics objects: muons (/i), 
electrons (e), photons (7), particle flow jets (anti-fcx)? 
and missing transverse energy (E™ 1SS ). 

Each event is sorted into precisely one event class, 
which represent a single final state, depending on its 
content of reconstructed objects. Together with the 



requirement of at least one charged electron or muon 
in any analysed event, this leads to about 250 event 
classes containing at least one event in data or simu- 
lation. 

Three kinematic distributions are examined, which 
are promising to spot new physics: 

• Scalar sum of the transverse momentum of all 
participating objects: ^Pt- 

• Combined (transverse) invariant mass M(x) of 
the objects in the event class. 

• E™ 1SS in classes containing E™ 1SS above our pre- 
defined threshold. 

Out of these distributions ^2pt is the most general 
observable, sensitive to many new physics models in- 
volving heavy new particles or modified high-energy 
behaviour. The invariant mass allows the discovery 
of new resonances. New physics with heavy or highly 
boosted "invisible" particles will show up in the E™ 1SS 
distribution. The kinematic variables are calculated 
from the objects passing the selection criteria listed 
below. Objects not fulfilling those criteria are not 
considered. 



2.1. Event and Object Selection 

Events are selected by single lepton triggers (muon 
or electron). The triggers require a minimal pt and 
additional quality criteria. Due to the rapidly chang- 
ing instantaneous luminosity, the thresholds have been 
increased over time, the highest being 15 GeV for 
muons and 22 GeV for electrons. In order to re- 
main well above the trigger thresholds, at least one 
muon with 25 GeV or one electron with 30 GeV is re- 
quired. The trigger efficiency on selected leptons is 
above 95%. 

A summary of the selection criteria applied on 
physics objects in each event is given in tab. [I[ 
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Table I: Selection criteria applied to the physics objects 
in each event. 



Object 


Pt 


M 


Other Cuts 


M- 


18GeV 


< 2.1 


isolation/track 


e 


25 GeV 


< 2.5 


shape /track / isolation 


y 


25 GeV 


< 1.4 


shape/isolation 


jet 


50 GeV 


< 2.5 


energy fraction 


7-imiss 


30 GeV 







To avoid using the same energy entry more than 
once, objects are removed if they are too close to each 
other (AR < 0.2) 1 : Jets are removed if there are pho- 
tons or electrons nearby, and photons are removed if 
an electron is close. 



3. Search Algorithm 



All kinematic distributions are fed through a scan- 
ning algorithm that systematically scans for devia- 
tions, comparing the simulated SM prediction with 
the measured data. Since all distributions are ana- 
lysed as binned histograms, the bin width is adjusted 
to the resolution of the considered variable. 

For each region of adjacent bins the algorithm cal- 
culates the number of expected events (£?), its uncer- 
tainty (<j), and the number of observed events (TV). 
The uncertainty calculation takes into account pos- 
sible correlations between bins and individual uncer- 
tainties (see section [5] for more details). A Poisson 
tail probability p can be computed which denotes the 
probability of a random fluctuation to be at least as 
extreme as the observed value. To incorporate the sys- 
tematic uncertainties on this probability, the Poisson 
distribution is convoluted with a normal distribution: 



p = < 



oo 

A - E 



-o-b) 2 e 



dfi if N < B 



(1) 



i=N 



e 2a 



d/i if N > B 



The factor A ensures normalisation, since the normal 
distribution is truncated at 0. In each distribution, 
the region with the smallest value of p = Pdata is cho- 
sen as the Region of Interest (Rol). 



4. Look-Elsewhere-Effect 



While Pdata denotes the probability of each Rol seen 
individually, it cannot be used as a statistical estima- 
tor for the global significance for such a deviation in 
any region. A penalty factor needs to be included to 
account for the number of investigated regions. Using 
(D(10 5 ) pseudo-experiments, we can compute a new 
probability p of measuring at least one deviation in 
any region in a given distribution with a lower p- value 
than the one seen in the Rol of the data. Each pseudo- 
experiment consists of one pseudo-data histogram for 
each distribution, thus representing one possible out- 
come of a measurement. The pseudo-data distribu- 
tion for each class is fed into the same scanning al- 
gorithm as described above, resulting in a collection 
of p- values p pse udo- The value p for a given distribu- 
tion is then simply the number of pseudo-experiments 
with Pp S eudo < Pdata, divided by the total number of 
pseudo-experiments : 



-^pseudo (.Ppseudo ^ Pdata) 



pseudo 



(2) 



5. Systematic Uncertainties 



This kind of search is highly affected by systematic 
uncertainties. An overview of the uncertainties used in 
this analysis is shown in tab. [TTJ Some of the included 
uncertainties can be correlated between bins of one 
distribution or across classes. 

Table II: Overview of the systematic uncertainties used 
in this analysis. 



Contribution 


Value 


Remarks 


MC statistics 


various 


sample dependent 


Luminosity 


4% 




parton density fkt. 


various 


PDF reweighting method 


jet energy corr. 


3 to 5 % 


p T and n dependent 


reconstruction eff. 


1 to 4 % 


object dependent 


misreconstrction prob. 


30 to 100 % 


object dependent 


W-boson cross sec. 


5% 


NNLO 


Drell-Yan cross sec. 


5% 


NNLO 


tt cross sec. 


10% 


NNLL 


di-boson cross sec. 


10% 


LO/NLO 


T cross sec. 


30% 


measured cross sections 


QCD-multijet cross sec. 


50% 


LO 


photon+jets cross sec. 


50% 


LO 



1 AR= y/Acf) 2 + At? 2 



The reconstruction efficiencies and the misidentifl- 
cation probability are determined from simulation and 
uncertainties are applied to cover possible differences 
to data. 
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6. Results 

As a demonstration two typical distribution are 
shown. The first one shows a Drell-Yan dominated 
class with two muons in the final state (fig. [I]) the sec- 
ond one shows a tt dominated class with one electron, 
one muon, two jets and missing transverse energy in 
the final state (fig. [2]). 

Overall 287 distribution in 118 event classes have 
been analysed in 2010 data. 
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Figure 1: A typical distribution as seen by the MUSiC 
analysis. The region in between the two blue dashed lines 
is the Region of Interest. The kinematic variable 
investigated here is the invariant mass of the two muons. 
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Figure 2: A typical tt dominated event class with one 
electron, one muon two jets and missing transverse 
energy in the final state. The kinematic variable 
investigated here is the scalar sum of the transverse 
momentum of all particles in the final state. 

The p distributions of the analysed event classes 
are determined separately for the three different kine- 
matic variables. 

As can be seen in fig.[3j|4j[5]the predictions for the p 
distributions are very accurate. The used SM samples 
describe the data very well. 

In the scope of this analysis no deviations from the 
SM have been found in LHC 2010 data. 
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Figure 3: Distribution of expected (light blue) and 
observed (crosses) p values for 2010 data for ^2pt- 
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Figure 4: Distribution of expected (light blue) and 
observed (crosses) p values for 2010 data for M( T ). 
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Figure 5: Distribution of expected (light blue) and 
observed (crosses) p values for 2010 data for E^ 1SS . 
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